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Abstract 

It is shown that the series derived by Mizrahi, giving the Husimi 
transform (or covariant symbol) of an operator product, is absolutely 
convergent for a large class of operators. In particular, the generalized 
Liouville equation, describing the time evolution of the Husimi func- 
tion, is absolutely convergent for a large class of Hamiltonians. By 
contrast, the series derived by Groenewold, giving the Weyl transform 
of an operator product, is often only asymptotic, or even undefined. 
The result is used to derive an alternative way of expressing expec- 
tation values in terms of the Husimi function. The advantage of this 
formula is that it applies in many of the cases where the anti-Husimi 
transform (or contravariant symbol) is so highly singular that it fails 
to exist as a tempered distribution. 
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1. Introduction 

A particularly useful and illuminating way of studying the classical limit is to 
formulate quantum mechanics in terms of phase space distributions 0, |[ ^. The 
advantage of such a formulation as compared with the standard Hilbert space for- 
mulation is that it puts quantum mechanics into a form which is similar to the 
probabilistic phase space formulation of classical mechanics. At least from a for- 
mal, mathematical point of view it thus allows one to regard quantum mechanics 
as a kind of generalized version of classical mechanics. 

There are, of course, many different phase space formulations of quantum me- 
chanics. The one which was discovered first is the formulation based on the Wigner 
function |, |, |. In the case of a system having one degree of freedom with 
position X, momentum p and density matrix p the Wigner function is defined by 

Wix,p)^^ J dye'Py{x-^\p\x + ^) 

(in units chosen such that h = 1). The Wigner function continues to find many 
important applications (in quantum tomography j^], for example). However, if 
the aim is specifically to represent quantum mechanics in a manner which resem- 
bles classical mechanics as closely as possible, then the Wigner function suffers 
from the serious disadvantage that it is not strictly non-negative (except in spe- 
cial cases 1^) — which makes the analogy with the classical phase space probability 
distribution somewhat strained. 

There has accordingly been some interest in the problem of constructing alter- 
native distributions, which are strictly non-negative, and which can be interpreted 
as probability density functions. There are, in fact, infinitely many such func- 
tions H, 1^, H, ^ 11 . The one which was discovered first, and which is the focus 
of this paper, is the Husimi, or Q-function §, 1 1, which is obtained 

from the Wigner function by smearing it with a Gaussian convolution: 

Q(a;,p) = i J dx'dp' exp [-{x - x')^ - (p - p'f]W{x' ,p') 

(in units such that h — 1, and where we assume that x, p have been made di- 
mensionless by choosing a suitable length scale A, and making the replacements 
x — > x/A, p — > Xp). It should be emphasised that it is not simply that the Husimi 
function has the mathematical significance of a probability density function. It 
also has this significance physically. It has been shown that the Husimi function 
is the probability distribution describing the outcome of a joint measurement of 
position and momentum in a number of particular cases |3[ ^, 18 . More 



generally it can be shown [|19|, ^ that the Husimi function has a universal signifi- 
cance: namely, it is the probability density function describing the outcome of any 



retrodictively optimal joint measurement process. In Appleby |20 it is argued that 
this means that the Husimi function may be regarded as the canonical quantum 
mechanical phase space probability distribution, which plays the same role in rela- 
tion to joint measurements of x and p as does the function \ {x \ "0)1^ in relation to 
single measurements of x only. 

If one wants to construct a systematic procedure for investigating the transition 
from quantal to classical it is not enough simply to find an analogue for the classical 
phase space probability distribution. One also needs an analogue of the classical 
Liouville equation, giving the time evolution of the probability distribution. In 
the formulation based on the Husimi function this is accomplished by means of 



Mizrahi's formula Il5 
see Lee [|, Cohen |j 



giving the Husimi transform of an operator product (also 
Prugovecki |||| and O'ConncU and Wigner 
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Let Aw denote the Weyl transform of the operator A, defined by ^, Q 

Awix,p) = J dye'Py{x-^\A\x+^) 

The Husimi transform (or covariant symbol) is then given by 

A^{x,p) = ^ J dxdp' exp [-{x - x)^ - [p-pf] Aw{x',p') (1) 

Mizrahi has derived the following formula for the Husimi transform of the 
product of two operators A, B: 

{AB)^ ^ A^e^+^-B^ (2) 

where d± = 2"^/^ (9a; =F idp). Using this formula, and the fact that the Husimi 
function is just the Husimi transform of the density matrix scaled by a factor 
l/(27r), it is straightforward to derive the following generalization of the Liouville 
equation: 

§iQ^{Hu,Q}n (3) 

where is the Husimi transform of the Hamiltonian, and {H^, Q}jj is the gener- 
alized Poisson bracket 

°° 2 

The first term in the sum on the right-hand side is just the ordinary Poisson bracket. 
The remaining terms represent quantum mechanical corrections. 

It is not apparent from Mizrahi's derivation, whether these expressions are exact, 
or whether they are only asymptotic. In Section S we will show that there is a large 
class of operators for which the series in Eq. (pf) [and consequently the series in 
Eq. (^J is absolutely convergent . This property is closely connected with the 
complex analytic properties of the Husimi transform, as discussed by Mehta and 
Sudarshan |^ and Appleby [ p6[ . 

The significance of the result proved in Section |2| is best appreciated if one 
compares Eq s. (2|4l|) with the corresponding formulae in the Wigner-Weyl formal- 



ism 




{AB)w = Aw exp | [^d^dp - dpdx) (5) 

d_ 

di' 



-W^{Hw,W}^ (6) 



where is the Weyl transform of the Hamiltonian, and where {i?Wj ^Iw denotes 
the Moyal bracket |2|. 

It can be seen that Eqs (||-^ and Eqs. (||-0) are formally very similar. However, 
this formal resemblance is somewhat deceptive, for it turns out that the two sets 
of equations have quite different convergence properties. 

The formula for the Weyl transform of an operator product, Eq. (^, is exact if 
either ^4 or i? is a polynomial in x and p (in which case the series terminates after 
a finite number of terms). More generally, if Aw, i?w are C°° functions satisfying 
appropriate conditions on their growth at infinity, then it can be shown that the 
series is asymptotic |30[ . However, there are many operators of physical interest for 
which the Weyl transform is only defined in a distributional sense, and for operators 
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such as this the series can be highly singular. Consider, for example, the parity 
operator V, whose action in the x-representation is given by 

(a:;| F I?/;) ~ (— x \tp') 

We have 

Vw{x,p) = 7t5{x)S{p) 
Substituting this expression into Eqs. (^) gives 

(V^)w{x,p) = 7rH{x)6{p)exp [| (S^^ - S^a^)] 6ix)5{p) 

The left-hand side of this equation — 1; whereas the expression on the right-hand 
side is an infinite sum each individual term of which is ill-defined (being a product 
of distributions concentrated at the origin) . 

Mizrahi's |lq] derivation of the formula for the Husimi transform of an operator 
product, Eq. (g), depends on the same kind of formal manipulation that is used 
in Groenewold's derivation of Eq. and so it might be supposed that the 
validity of the formula is similarly restricted. However, it turns out that the sum 
in Eq. (||) is actually much better behaved. In fact, it will be shown in Section ^ 
that, subject to certain not very restrictive conditions on the operators A and 13, 
the sum on the right-hand side of Eq. (||) is not only defined and asymptotic; it is 
even absolutely convergent for all x, p. This is essentially because is typically 
a much less singular object than [due to the Gaussian convolution in Eq. (^]. 

In Section ^ we apply the result just described to the problem of expressing 
expectation values in terms of the Husimi function. 

The expectation value of an operator A can be obtained from the Wigner function 
using the formula ||l|, ^ 

Tr(/3i) = J dxdpAw{x,p)W{x,p) (8) 

In certain cases we can also express the expectation value in terms of the Husimi 
function using ||l|, 

Tt{pA) = J dxdpA^{x,p)Q{x,p) (9) 

where A-^ is the anti-Husimi transform (or contravariant symbol) of A, defined by 

(10) 

[with d± ~ 2^^/'^{dx ^idp), as before]. Eq. (^) is valid (for example) whenever ]3l|| 
exists as a tempered distribution and Q belongs to the corresponding space 
of test functions {i.e. the C°° functions of rapid decrease). However, we have the 
problem that A^ is often so highly singular that it is not defined as a tempered 
distribution — which means that the usefulness of Eq. (^ is somewhat limited. This 
is often seen as a serious drawback of the Husimi formalism. 

However, it turns out that it is often possible to circumvent this difficulty. Sup- 
pose we substitute the series given by Eq. ( p^ ) into the right hand side of Eq. (||) , 
and suppose we then reverse the order of sum and integral. This gives 

Tr(M) = E / "^""^^ {dld-A^{^,P)) Qi^^P) (H) 

In Section ^ we show that it often happens that the sum on the right-hand side of 
this equation is absolutely convergent, even in many of the cases where Ajj- fails to 
exist as a tempered distribution. 
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2. Convergence of the Product Formula 
We will find it convenient to work in terms of coherent states. Define 

and let </>„ denote the n^^ (normalised) eigenstate of the number operator 
00 = (/)„ = ^(a'l')"0o 



a' a: 



"n 

Let Dxp be the displacement operator 

n _ i{px-xp) 

and define 

The (j)xp are the coherent states. Let A be any operator (not necessarily bounded) 
with domain of definition D^, and suppose that 4>xp € for all x, p. It is then 
straightforward to show that 

^iiix,p) = {<l^xp, A(l)xp) 

(we no longer use the Dirac bra-ket notation, because the existence of {ip, Ax) does 
not, in general, imply the existence of {A'^ip, x))- 

If A, B are both bounded then the proof of Eq. (H) is comparatively straightfor- 
ward. However, we want to make the proof as general as possible. We then have 
the difficulty that the sum in Eq. (^) will only be defined if Ajj and _Bjj are both 
C°°] whereas functions of the form {Dxpip, AD^pip) are, in general, not even once- 
differentiable, let alone C°°. We are thus faced with the question: what conditions 
must we impose on the operator A in order to ensure that the function A^ is C°°7 
One answer to this question is given by the following theorem. 

Theorem 1. Let 2)^, T>^^ be the domains of definition of A, A^ respectively. 
Suppose that (pxp G 2)^ H D^t for all x, p. Suppose, also, that {(pxipi, A(j)x2P2)j 
{4'i;xipn A(j)x2P2) ^'^'^ {4>i;xipn -A^ 4>X2P2) ^'"^ continuous functions onW^. Then A^ 
is an analytic function, which uniquely continues to a holomorphic function defined 
on the whole of . 

The continuation is given by 

A^{x,p) = (^-P-'-^^-.P.) (12) 

where x,p are arbitrary complex, and where x±,p± are the real variables defined by 

x± = ^{x + x*)±'-{p~P*) (13) 

P± = l{p + P*)T^{x-x*) (14) 

This theorem is a strengthened version of results proved by Mehta and Sudar- 
shan and Appleby ||2^ . The proof is given in Appendix 

It is worth noting that the condition in the statement of this theorem is quite 
weak. If the three functions listed exist and are continuous then, without making 
any explicit assumption regarding the differentiability of these functions, it auto- 
matically follows that Ajj must be complex analytic. 

We also have the following lemma: 
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Lemma 2. Suppose that A satisfies the conditions of Theorem Th 

QT 



{<t>n-x. 
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1 
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n 



(^+-^-)""'— ^h(^,p) (15) 



-4)"-''^^h(^,p) (16) 



r=0 

where x±, p± are the variables defined by Eqs. ([I^j and (jZ^, and where 

z± = :^(^=^*p) 

The proof of this lemma is given in Appendi x |B| . 

If X, p are both real (so that Z- — z*^_) Eqs. ( |l5|) and ([l6| ) become 

'-^n^xp, i-Z-xp) = ^ 5!!Ah(x,p) (17) 



where 



(i^xp, <^„;xp) = 4=f ^?^H(a;-P) (18) 
n 



d - — - — f— I— 

dz± y/2 \dx dp 

Using these results the proof of the product formula becomes very straightfor- 
ward. Let A, B be any pair of operators satisfying the conditions of Theorem ^. 
Suppose, also, that <j)xj} G ab '^'^'^ '"^^^ G R. Then, for real x, p, 

{Ab)^{x,p) = {(f>xp, AI3(f>xp) 

= {A^(t>xp, B4>xp) 

OG 

) (19) 

where we have used the fact that the (j)n;xp constitute an orthonormal basis. The 
absolute convergence of this sum is an immediate consequence of basic Hilbert space 
theory. Using Eqs. (^^ and ( [l8| ) we deduce 

{AB)^{x,p) = ^dlA^{x,p)d^B^{x,p) = A^{x,p) e^^B^{x,p) 

which is the product formula. The absolute convergence of this sum follows from 
the absolute convergence of the sum in Eq. (|l9|). 

3. Expectation Values 

We now discuss the implications that the result just proved has for the con- 
vergence of Eq. (pi]), giving the expectation value of A in terms of the Husimi 
function. 

Of course, one does not expect the right hand side of Eq. (|l^) to converge for 
arbitrary A and p — since, apart from anything else, an unbounded operator does 
not have a well-defined expectation value for every state p. We therefore need to 
place some kind of restriction on the class of operators A and density matrices p 
considered. The result we prove is probably not the most general possible. However, 
it will serve to illustrate the point, that the sum on the right hand side of Eq. ( pT| ) 
is often absolutely convergent, even in many of the cases where the anti-Husimi 
transform fails to exist as a tempered distribution. 
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We accordingly confine ourselves to the case of density matrices for which the 
Husimi function Q e J(Ili^), where II(M^) is the space of C°° functions which are 
rapidly decreasing at infinity [Slj l {i.e. the space of test functions for the space of 
tempered distributions). In other words, we assume that 

sup \il + x^ +pydl''d;Qix,p)\<oo 

(3:,p)eR2 

for every triplet of non-negative integers I, m, n. 
We assume that A has the properties 

1. (t>a;p e Tiji^r] 2)^-1. for aU x,p G K. 

2. There exist positive constants K± and non- negative integers N± such that 

\\A4>,p\\<K^{l+x^+pY- (20) 
\\AU.p\\<K+{1 + x^+pX^ (21) 

for all x, p e M. 

We will say that an operator satisfying these two conditions is polynomial bounded. 
The following lemma gives two properties of such operators which will be needed 
in the sequel. 

Lemma 3. Suppose that A is polynomial bounded. Then 

1. A satisfies the conditions of Theorem^. In particular, A^ is analytic. 

2. For every pair of non-negative integers m, n there exists a positive constant 
Kmn o,nd a non-negative integer Nmn such that 

\d:'d;A^{x,p)\ < K^,,{l + x^+p^)''- (22) 

for all a;, p e M. 

The proof is given in Appendix |^. 

We now ready to prove the main result of this section. 

Theorem 4. Suppose that A is polynomially bounded, and suppose that the density 
matrix p is such that the corresponding Husimi function is rapidly decreasing at 
infinity. Suppose, also, that Ap is of trace-class. Then 

Tr(ip) = ^ ^—J- / dxdp {dld^A^{x,p)) Q{x,p) (23) 

where the sum on the right hand side is absolutely convergent. 
Proof. We have 

Tr(i/5) = ^ / dxdp {(j>xp, Ap(f)^p) 



= ^ J dxdp ( ^(^Va;p, 4'n;xp) {4>7i: 



Xp7 f-"VXp/ 



We now use Lebesgue's dominated convergence theorem to show that we may 
reverse the order of sum and integral. In fact, it follows from the Schwartz inequality 
that 



^^(^A^ (pxp^ 'Pn-jXp) {'Pn;xp^ p4'xp) 



n=0 



< 



< II^Vxpll WP'f'xpl 



El 

n=0 
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We have 

WpipxpW = {{<Pxp, P'^<l>xp))^ < {{(j)xp, P(j)xp)Y = yj2-KQ{x,p) 



which, together with Inequality (21), implies 



By assumption, Q{x,p) E !}(M^). It follows that ||ylt(/)2.p|| ||p(/)^p|| is integrable. We 
may therefore use Lebesgue's dominated convergence theorem to deduce 

1 °° /■ 

Tr(i/5) = ^ X! / dxdp {A^ (j>j;p, (j>n-xp){(l>n;xp, P(l>xp) (24) 



n=0 

where the sum is absolutely convergent, since 



dxdp {A^ (fixp, 4>n;xp){4'n;xp, HVxpi 



< I dxdpWA"^ 4>xp\\ \\p4>xp\\ < 



E 

ri=0 

We know from Lemma § that A satisfies the conditions of Theorem |l|. We may 
therefore use the results proved in the last section to rewrite Eq. ( ^4| ) in the form 

Tr(ip) = ^l f dxdpdlA^{x,p)dlQ{x,p) 

n=0 ' 

Finally, it follows from Inequality (|2^), together with the fact that Q E 3(R^), that 
we may partially integrate term-by-term to obtain 

Tr(ip) = J2 —T- / dxdp{dld^A^ix,p))Qix,p) 



□ 



The right hand side of Eq. rtq) (expressing (A) in terms of the Wigner function) 
is defined whenever W E 3{R^ and Ayj exists as a tempered distribution. On 
the other hand, although it is true that Q E whenever W E (see 

Theorem IX. 3 of ref. the fact that exists as a tempered distribution is not 

evidently sufficient to ensure that A is polynomial bounded. So we have not shown 
that Eq. (|23| ) has the same range of validity as Eq. However, it can be shown 
that A is polynomially bounded if <j)xp E n D^t, and if {A''A)^ff and {AA^)^^ 
exist as tempered distributions (see Theorem IX. 4 of ref . [0). In the applications 
one meets with operators satisfying these conditions much more commonly than one 
meets with operators for which Aj^ exists as a tempered distribution. For instance, 
every bounded operator is polynomially bounded, whereas there are many bounded 
operators of physical interest for which A^ fails to exist as a tempered distribution. 
The above result consequently represents a significant improvement on the results 
that were previously known. 

Of course, just from the fact that Eq. (|2^ ) is convergent, it does not necessar- 
ily follow that the convergence is sufficiently rapid to make the formula useful in 
practical, numerical work. This question requires further investigation. 

4. Conclusion 

As has been stressed by Mizrahi , Lalovic et al |^ , Davidovic and Lalovic 
and others, the Husimi formalism provides an especially perspicuous method for 
studying the relationship between quantum and classical mechanics. It establishes 
a one-to-one correspondence between the basic equations of the two theories, so 
that one can start with a classical formula, and then turn it into the corresponding 
quantum formula by adding successive correction terms. Moreover, the fact that 
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Q{x,p) describes the outcome of a retrodictively optimal joint measm'ement of x and 
p means that one could reasonably argue that the Husimi function is the 

most natural choice for a quantum mechanical analogue of the classical probability 
distribution. 

In this paper we have investigated the convergence properties of two of the key 
formulae in the Husimi formalism. We have shown that the formula giving the 
Husimi transform of an operator product has much better convergence properties 
than the corresponding formula in the Wigner function formalism. In particular, 
the Husimi formalism leads to a convergent generalization of the Liouville equation 
for a very large class of Hamiltonians. We have also shown that the convergence 
properties of the formula expressing the expectation value (A) in terms of the 
Husimi function, although seemingly not as good as those of the corresponding 
formula in the Wigner function formalism, are significantly better than the often 
highly singular character of Ajj would suggest. 

These results lend additional support to the suggestion that, in so far as the 
aim is specifically to formulate quantum mechanics as a kind of generalized version 
of classical mechanics, then the formalism based on the Husimi function has some 
significant advantages. 

Appendix A. Proof of Theorem |i] 
For arbitrary complex x, p define 



Fix,p) 



where x±, p± are the (real) variables defined by Eqs. ( p^ ) and (14). 
It is easily seen that, if x, p are both real, then 

F{x,p) = {(j)xp, A(j)xp) = A^{x,p) 

The problem thus reduces to that of showing that, if A has the properties stipulated, 
then F is holomorphic. We will do this by showing that F satisfies the Cauchy- 
Riemann equations with respect to the complex variables 

z± ^ -^{x± ± ip±) ^ ^{x ± ip) (25) 

In fact, it is straightforward to show that (j)xp, regarded as a vector- valued func- 
tion of two real variables, is differentiable in the norm topology; the derivatives 
being given by 



d . 1 



dx ~ ^J2 2^^^P 
9 _ i i 
'dp ~ 2^'' 



xp 



Also, 



exp[-^(a;+ - x^f - ^{p+ - p^f + i^{p+X- - P-x+)\ 



Consequently, F is differentiable with respect to the variables X-, p-. Moreover 

dx^ ' ~ 72 (0 ' V2^ - V 

^i-^F{x,p) (26) 

from which it follows that F satisfies the Cauchy-Riemann equations with respect 
to the complex variable z„ = (a;_ — ipJ)j \f2. 
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We can alternatively write 

Consequently, F is also difFerentiable with respect to the real variables a;+, p+. 
Moreover 



dx+ ' V2 ) V2 

= -i-^F{x,p) (27) 

from which it follows that F satisfies the Cauchy-Riemann equations with respect 
to the complex variable z+ = (a;+ + ip+)/^/2. 

If A has the properties specified in the statement of theorem, then we see from 
Eqs. ( p6|) and ( |27| ) that the partial derivatives dF/dx±, dF/dp±, are continuous 
functions on ]R\It follows that F is a holomorphic function of the complex 
variables z±. Referring to Eq. ( p5| ) it can be seen that the variables z± are linear 
combinations of x, p. We conclude that is a holomorphic function of x, p. 

Appendix B. Proof of Lemma H] 

It is straightforward to show that (f>n;x^p- , regarded as a vector valued function 
of two real variables, is differentiable in the norm topology. Moreover 

1/9 . a \ , , 1 



Hence 



y/2 \dx- dp- 



A4>x+p^) / ^ \ (<^n;x_p_, 



+ 1 7- T + {Z_ - Z+) — ^ 

\(Px-P-j Vx^p^l \(Px-P-y (Px^p^l 



Iterating this result, and using 



d d Y d"^ 

we obtain Eq. (|l5|). 

The proof of Eq. ( [l6| ) is similar. 

Appendix C. Proof of Lemma |3| 

Proof of (1). We need to show that, if A is polynomial bounded, then the functions 
{4>xipt, A(j)^^p^), (0i;xipi, A(j)x2P2) and {(t)i-xipt, i-Vajaps) are continuous. 
Consider the function {(j^i-^xipn ^4'x2P2)- We have 

^l\x'^p^^: '^4'x'2P'2) {01;a;ipii ^0a:;2P2 ) 

< (('/>l;i;'ip'i - ^l;2:ipi), ^(Axip'J + (01;xipi , ^((/>2:'2Pi ^ ^2:2P2 )) 

In view of Eq. (|20| ) we have 

((01;2:'iP'i - 01;xipi), Mx'2P'2) < K^{l + x'2 + ' \\4>Ux'^pi " 01;2:ipill 
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Also, using the completeness relation for coherent states, together with Eq. (|21|), 
we find 



1 

< f{xi,Pl)\\(t)x'2P'2 - (t'x2P2\\ 

where / is the polynomial 
K f 



2i7r 



dx'^dp's ^ I + x'i + p':,^l + (4 + xif + [p'^+Pif) 



2\N+ 



X exp 



Putting these results together we find 

< K^{1 + X2 + P2)^'\\(t>l;x'^p'^ - (t)l;xipA\ + /(xi , pi ) || " 0a;2P2 II 



^xp 



and 4'i;xp, regarded as vector- valued functions on R'^, are continuous in the 
norm topology. Consequently 

(01;x'ip'i, Mx'2P'^) {(f>l;xipi, Mx2P2) 

as {x[,p'i, x'2,P2) {^iiPii ^2tP2}- It follows that {ipi-xipi, ^0X2P2) continuous. 
Continuity of the functions {(f^xipi, -^4'x2P2} ^^'^ {<l>i;xipi, 4'x2P2) proved in the 
same way. 

Proof of (2). The completeness relation for coherent states implies 

1 



27r 



dx'dp' {4>xp, 4'x'p'){A^4'x'p', 4>xp) 



Using 



d+(l)n-xp = Vn+l(j>^n+l);xp + -ZZ* 



9_ 



2 ^4^xp 



if n = 
if n > 



1 \/n(f'{n — l);xp 2z4'n;xp 

[where d± — 2^^/'^{dx T idp) and z — 2^^/'^{x + ip))], and differentiating under the 
integral sign, it is not difficult to show that 

n+m „ 

d'^dpA^{x,p) = Crs dx'dp {4>r:xp, <t>x' p') {A^ (px' p' , <l>s;xp) 

r,s=0 •' 

for suitable constants Crs- We have 
1 



22 Vr-! 



((x'-x)2 + (p'-p)2)^cxp 



-h^x -xf - ^{p -pf 



In view of Inequality ( pT| ) it follows that 

/ dx'dp' {(pr-xp, (t)x'p'){A^4'x'p', (f>s;xp) 
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< 



I 



dx"dp"{l + {x" + xf + {p"+pf) 



X [x" +p" j^exp -^[x +P ) 



It can be seen that the expression on the right hand side of this inequahty is a 
polynomial in x and p. Consequently 
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